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In this paper we consider dynamics of a molecular system subjected to external pumping by a 
light source. Within a completely quantum mechanical treatment, we derive a general formula, 
which enables to asses effects of different light properties on the photo- induced dynamics of a molec- 
ular system. We show that once the properties of light are known in terms of certain two-point 
correlation function, the only information needed to reconstruct the system dynamics is the reduced 
evolution superoperator. The later quantity is in principle accessible through ultrafast non-linear 
spectroscopy. Considering a direct excitation of a small molecular antenna by incoherent light we 
find that excitation of coherences is possible due to overlap of homogeneous line shapes associated 
with different excitonic states. In Markov and secular approximations, the amount of coherence 
is significant only under fast relaxation, and both the populations and coherences between exciton 
states become static at long time. We also study the case when the excitation of a photosynthetic 
complex is mediated by a mesoscopic system. We find that such case can be treated by the same 
formalism with a special correlation function characterizing ultrafast fiuctuations of the mesoscopic 
system. We discuss bacterial chlorosom as an example of such a mesoscopic mediator and propose 
that the properties of energy transferring chromophore-protein complexes might be specially tuned 
for the fluctuation properties of their associated antennae. 



I. INTRODUCTION 

In recent years, primary processes in photosynthesis 
have received a renewed interest from a broader physi- 
cal community thanks to experimental observation of co- 
herent energy transfer in some photosynthetic systems. 
The ground breaking coherent two-dimensional electronic 
spectroscopy (2D-ES) experiment of Engel et al. [l| has 
led to new appreciations of the role that may be played by 
coherent dynamics in excitation energy transfer (EET), 
and of the quantum mechanical nature of photosynthetic 
systems in general [2']. Special theoretical effort has been 
made to understand the role of noise 0-0| in the dynam- 
ics of excitation energy transfer, and the role of coher- 
ence jM^Ill ill excitonicaly coupled systems. On the ex- 
perimental front, the method of coherent 2D-ES fl3l, [T^ 
has established itself as a tool opening new window into 
the details of energy transfer dynamics in photosynthetic 
[T5l - [20| . and other molecular systems 0, |2l|,[22|. Coher- 
ent effects have been now reported in different molecular 
systems, often biologically relevant [l^, - a general- 
ity that asks for a search of the possible evolutionary 
advantage underlying their abundance in photosynthetic 
pigment-protein complexes. 

The principle pigment molecules responsible for the 
primary processes of photosynthesis are chlorophylls 
(Chls) and bacteriochlorophylls (BChls) [13, [11. They 
are involved in accumulation of light energy via the exci- 
tation energy transfer to specific pigment-protein com- 
plexes - reaction centers. Spectral variability of pho- 
tosynthetic light-harvesting pigment-protein complexes 
arises either from excitonic interactions between pigment 



molecules or from their interactions with protein sur- 
rounding. Both these interactions are the main factors 
determining the excitation dynamics in light-harvesting 
[2^. Excitonic aggregates are subject to interaction with 
two types of environments, and they provide means of 
transferring energy from one environment to another. 
First of these environments, the radiation, is under nat- 
ural conditions at much higher temperature than the 
second environment, the protein scaffold and indeed the 
photosynthetic chemical machinery as a whole. The ex- 
cess of photons on suitable wavelength in the radiational 
environment is used to excite spatially extended antenna 
systems that concentrate excitation energy to the reac- 
tion center, which in turn drives charge transfer processes 
across cellular membranes to create the transmembrane 
potential and the pH gradient [2^ . 

Non-equilibrium processes occurring in photosynthetic 
systems during light harvesting are conveniently de- 
scribed by reduced density matrix (RDM) theory (26l - 
[28| which has an advantage of being applicable to dis- 
ordered statistical ensembles that the experiments often 
deal with. However, with recent 2D experiments that 
enable us to distinguish the homogeneous and inhomoge- 
neous spectral broadening, and with the progress in single 
molecular spectroscopy [2^ we can gain insight into the 
time evolution characteristic to single molecules interact- 
ing with their environment (sol . [si ] . This fact enables us 
to return to the wavefunction formalism and to look at 
light harvesting from the point of view which takes the 
superposition principle of quantum mechanics seriously. 
It has been show that such an approach yields many in- 
teresting insights into the emergence of the classical prop- 
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erties of molecular system from their underlying quantum 
mechanical nature |32l. |33|. As the light-harvesting pro- 
cesses seem to operate on the interface between classical 
and quantum worlds it seems appropriate to look at them 
from the point of view of the decoherence program of Zeh, 
Zurek and others (s^. l35l|. 

The process of light harvesting could then be describes 
as follows. First, the system is in an "equilibrium" initial 
state l^'o) characterized by the excitonic ground state 
\g), the state of protein (phonon) environment |$p) cor- 
responding to this electronic ground state and some state 
of light |So), i.e. 

|vI/o)-|5>|$b)|So>. (1) 

The light-harvesting occurs when the state of light is such 
that the time evolution of the system leads to population 
of higher excited states |e„) of photosynthetic antenna. 
These states are formed from excited states of Chls and 
other chromophores, such as carotenoids [1^. We denote 
these combined excited states as excitons. In the first ap- 
proximation, photosynthetic antenna remains in the ex- 
cited state until the excitation energy is transferred to the 
reaction center. This happens much faster than compet- 
ing process of spontaneous emission which can therefore 
be neglected in our discussion. When the interaction of 
the antenna with light is switched on, the change occur- 
ring in the ground state portion of the total state vector 
after the passage of time At is 

I *o) ^ «At |5) |$B> |So) + ^ /3i") |e„) |<i>B) (2) 

n 

The subsequent time evolution of the excited state por- 
tion of the state vector is independent of the ground state 
part, and we can thus look at it separately. Because we 
neglect spontaneous emission, any excitation to higher 
excited state, as well as transitions between exciton states 
due to the light, the state vector remains approxi- 
mately uncntangled with excitons and the protein bath 
for the rest of the energy transfer process. It can there- 
fore be omitted. The initial state for the energy transfer 
process thus reads 

\^,ito))=J2l3^-Hto)K)\^B), (3) 

n 

where we omitted the lower index At. If the basis of the 
states |e„) is chosen so that the molecular Hamiltonian 
is diagonal, the energy transfer occurs only due to inter- 
action of excitons with their surrounding environment. 
This interaction leads to an entanglement of excitons and 
the environment 

|*,(t))=^/3(")(t)|e„)|<i>i,")(t)). (4) 

n 

After a sufficiently long time the environment state vec- 
tors corresponding to different electronic state diverge 



maximally and the reduced density matrix becomes di- 
agonal in some basis, i.e. 

p(i)=irs{|*e(<))(*eW|} 

= 5] (t) it))* (<i>i;"' (t) I (t)). (5) 

mn 

Often, to a good approximation, such preferred basis is 
the one in which the electronic Hamiltonian is diagonal, 
the so-called excitonic basis. However, notable correc- 
tions to this rule are predicted even for weak system-bath 
coupling [l^ [3^. 

The final state of the energy transfer is the one in which 
just reaction centers are populated 

l^^it)) ^/3(«'^'=)(t)|efic.)l4''''^'W)- (6) 

k 

The last step of the energy transfer, from the antenna 
to the reaction center is often slower than typical trans- 
fer times between antenna complexes, and so the final 
state is well localized on the reaction center, and coher- 
ences between individual reaction centers are unlikely to 
survive. 

It is clear from the above discussion, that decoherence 
during the energy transfer in the antenna is determined 
by the evolution of the environmental degrees of free- 
dom (DOF). The decoherence from the rest of the sys- 
tem might be required for the localization of the energy 
in the reaction center, but there is no obvious reason for 
fast decoherence during the initial steps of energy transfer 
in the antenna, apart from the fact that a bath formed 
by a completely random disordered environment would 
lead to just such fast decoherence. It has been suggested 
before that the protein environment might play a more 
active role in steering and protecting electronic excita- 
tion [11, [2^ and controlling the decoherence might be one 
of the possible pathways to more robust EET. 

There is however one important caveat in the above 
scheme. The initial condition, Eq. ([3]), has been intro- 
duced artificially into Eq. ([5]) as a result of an interac- 
tion occurring during some short time interval At. If the 
system is continuously pumped, individual contributions 
similar to Eq. ^ will interfere, possibly disabling any 
effect of cooperative involvement of the bath. It is even 
more important to consider the question what are the 
effects of natural sun light [13, i-e. whether the coher- 
ent scenario outlined above is plausible for the photosyn- 
thetic system in vivo, or not. This depends strongly on 
the nature of the excitation process, whether it occurs in 
discrete independent jumps of the kind described by Eq. 
(12) , or continuously over a long period of uncertainty in- 
terval of the photon arrival. The former view is usually 
held in support of the relevance of ultrafast spectroscopic 
finding for in vivo function of the photosynthetic systems 
[m . Below we derive a general formula which enables us 
to describe all these regimes by a unified formalism, and 
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also enables us to place the observables of ultrafast coher- 
ent spectroscopy in perspective with the dynamics under 
natural conditions. In a somewhat extended form our re- 
sult is also applied to another case cited in support of the 
utility of coherent dynamics in photosynthetic systems, 
a case where a small photosynthetic complex is excited 
through another, possibly mesoscopic, antenna (TT| . 

The paper is organized as follows. Next section in- 
troduces a rather general model of photosynthetic aggre- 
gate. In Section Hn] we discuss the dynamics of a system 
excited by coherent pulsed light and the observables of 
the ultra- fast no n- linear spectroscopy. Section ITVl is con- 
cerned with the excitation of a photosynthetic system by 
the light from a general source. Implications of the theory 
for excitation by thermal and coherent light, as well as 
excitation mediated by mesoscopic system are discussed 
m Section El 



II. HAMILTONIAN OF A MODEL 
PHOTOSYNTHETIC SYSTEM 

In this section we briefly review the excitonic model 
that was very successfully applied to model the spec- 
troscopic properties of Chi- and Bchl- based light har- 
vesting chromophore-protein complexes (see e.g. [l6j). 
We assume N monomers with ground states \gn), excited 
states |e„), n = 1, . . . , A^, and with electronic transition 
energies e„. These monomers are interacting with the 
phonon bath of protein DOF so that the Hamiltonian of 
the monomer reads 

H,, = {T + Vg)\gn){9n\ + (in +T+ Ve)\in){en\- (7) 

Here, T is the kinetic energy operator of the bath, and 
Vg and Ve are the potential energy operators of the bath 
when the system is in the electronic ground- and excited 
states, respectively. We set the ground state electronic 
energy to zero for conveniency. The Hamiltonian, Eq. 
([7]) can be split into the pure bath, pure electronic and 
the interaction terms so that 

Hn = {T + Vg) ® 111 + Is ® (en + {Ve " Vg) ,,)\er,) {in\ 

V ' v ' 

+ {Ve-Vg-{Ve-Vg),g)^\er,){er,\. (8) 
^ ^ ^ 

Here, Is is the unity operator on the bath Hilbert space 
and Im is the unity operator on the Hilbert space of the 
electronic states. The equlibrium average {Ve — Vg)eq 
of the potential energy operators was added to the elec- 
tronic energy so that the interaction term is zero for the 
system in equilibrium. 

In chromophore-protein complexes many such 
monomers are coupled by resonance coupling. The 



whole complex can be described by means of collective 
states including the ground state 

N 

\9)=X{®\9n), (9) 

n=l 

one excitation states 

a-l N 

|ea) = n ^Ifi'") ® ® n ®\9yn), (10) 

n—1 m—a-\-l 

and states containing higher number of excitations. For 
the sake of brevity we now stop writing the symbol of 
the direct product ® and the unity operators etc. ex- 
plicitely. The total Hamiltonian of the complex including 
resonance interaction is then defined as 

N n 

Hb + Hm + Hm-b = ^ + ^ H]^j 

n—1 n—1 
N n 

B- (11) 

n^m n—1 

If the system-bath interaction with bath is weak, the re- 
ferred basis into which the electronic system relaxes due 
to interaction with the bath is, to a good approximation, 
the one in which the electronic part of the Hamiltonian 
is diagonal. Let us denote these states as |e„). They are 
usually termed excitons and they represent certain linear 
combination of the collective states |e„) where excitations 
are localized on individual chromophore molecules. One 
of the most important characteristics of this model is that 
it does not include direct relaxation of the electronic ex- 
cited states to the ground states due to electron-phonon 
coupling. This is well satisfied by Chls and BChls on the 
ultrafast time scale of which light harvesting processes 
occur. 



III. EXCITATION BY COHERENT PULSED 
LIGHT AND NON-LINEAR SPECTROSCOPY 

Let us now consider experimental methods which pro- 
vide information about time evolution of excited states of 
photosynthetic systems. Because of the timescale of EET 
processes, spectroscopy with ultrashort time resolution is 
a necessary tool. The interaction of the pulsed coherent 
light with the photosynthetic system is well described in 
semi-classical approximation [3g. Electric field of the 
light is then considered as an external parameter of the 
system Hamiltonian. Electronic DOF can be prepared 
very fast in an excited state, not affecting, to a good ap- 
proximation, the bath DOF. Thus, in an experiment with 
an ideal time resolution, we would have the system pre- 
pared in the excited state, Eq. ([3]). The time evolution 
of the system is governed by the Schrodinger equation 

^14-^(0) = -\{Hb + Hm + i/M-B)|*e(t)) 
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+ 5(t)|*e(io)), (12) 

with initial condition \'i'e{t)) ~ for t < to. The last 
term in Eq. describes the ultrafast event of the 

molecule-radiation interaction. Formal solution of this 
equation reads |*e(i)) = C/s(i)C^M(i)t^M-s(OIV'e(io)), 
where we defined evolution operators C/_b(^), UM{t) of 
the bath and the molecule, respectively, as 

Ueit) = e{t-to)exp{-^HB{t-to)}, (13) 



Unit) = e{t - to) cxp{--HM{t - to)}, (14) 

and the remaining interaction evolution operator as 

t 

UM-sit) = Q{t-to)exp{-'-l dTUl(T)uliT) 

to 

X HM-BUMiT)UBiT)}, (15) 

After excitation, the process of energy transfer proceeds 
according to the description presented in Introduction 
and can be experimentally monitored. 

A. Evolution superoperator 

Matrix elements of the RDM of the molecule, which 
holds the information about the population probabilities 
and the amount of coherence between electronic states 
are given by expectation value of projectors |e„)(em|, 

Pnm 

(t) = {Mt)\e^){en\Mt)) 

= trB{{en\i'eit)){^eit)\ern)} 

= {en\trB{UM {t)UM-B it) J2 (Z?^'^ )* l«) {b\ 

ab 



X \^B)i^Ul_s{t)Ul{t)}\ern) (16) 

This can be rewritten by defining an evolution superop- 
erator U{t) which acts on initial density matrix poWeq, 
i.e. 

{t)^trB{{en\W{t)\e^)} 



= (e„|Z^('=)(i)PoWe,|e™), (17) 



The matrix elements of the superoperator read 

U^iS)^{<^\UMmM-Bit)\c)... 

X {d\ui^smUm, (18) 

where the dots . . . denote where an operator on which 
U'^'^\t) acts has to be inserted. The reduced evolution 
operator U''^'^ (t) defined as 

U^^\t)=trB{U^''\t)}, (19) 

contains information about the evolution of the RDM 
only. 

B. Non-linear spectroscopy 

In non-linear spectroscopy, coherent laser light is used 
to investigate the dynamics of molecular systems by ap- 
plying special sequences of pulses. Some pulses act to in- 
duce non-equilibrium dynamics (pump) , and other pulses 
act to monitor (probe) the evolution after the pump. One 
of the most advanced of these methods, coherent 2D-ES 
[13, [HI measures the response of a system to three pulses 
traveling in different directions fci , k2 and fcs . The detec- 
tion is arranged in such a way (measuring in the direction 
— fci + k2 + k^) that the signal is predominantly of the 
third order, with contributions of one order per pulse HI] . 
Let us denote delays between the first two pulses by r and 
the delay between the second and the third pulse by T. 
If the pulses are ideally short, the signal is composed of 
two kinds of contribution. First, contribution that in- 
volves population of the excited state corresponding to 
the density operator 

<k(*'^) = I^^^W)(V4'=^H^ + ^)l' (20) 

and second, contribution that involves evolution in the 
ground state 

witSt^r) = |*B)|g)(^f-^^Hi,t + r)|. (21) 

Here, we denote the pulses acting on the state vector 
by their corresponding wave vector in the upper index, 
and the excited state or ground state bands by the lower 
index g and e, respectively. For these statistical operators 
we can define evolution superoperators U'^'^\t), U'^^^t), 
U'^^3){^t) and [/(^^^(t) in analogy with Eqs. ^ and P^ . 
so that 

wi%Si,T)^U^^\t)U^^^\T)\^B){^B\ti\9){g\n, (22) 
and 

wl^l^^{t,T)^U^^\t)U^^^\T)\^B){^B\p'\9){9\. (23) 

The superoperator hl'^'^^^ (t) is the evolution superopera- 
tor of a coherence projector |en)(5| and analogically 
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for U^^'^\t). After a delay T the third ultrafast pulse is 
applied and the non-linear signal is recorded. The signal 
corresponds to indirectly to non-linear polarization of the 
sample, and is usually measured in frequency domain 

(24) 

Here, we denoted 

oo 

ZY("»)(w) =y"dt e*"*W('=s)(t), (25) 



and Pg — \g)(^g\. In 2D coherent spectroscopy, the signal 
is in addition Fourier transformed along the time delay 
T, so that the spectrum is defined as 



oo 




'-^i;;(3)(w,,r,r). (26) 



— oo 



The spectrum defined in this way has a suitable inter- 
pretation of an absorption - absorption and absorption 
- stimulated emission correlation plot, with a different 
waiting times T between the two events. The 2D spec- 
trum is in practice measured with finite pulses, and the 
measured time domain signal is thus a triple convolution 
of the responses to a delta pulse excitation, with the ac- 
tual finite pulses (4(]j . 

From this rough sketch of the principles and the infor- 
mation content of the coherent 2D spectroscopy it should 
be clear that 2D spectroscopy is aimed at disentangling 
the dynamics of the system during the time delay T . In 
the so-called Markov approximation, when the dynamics 
in time intervals t, T and t is assumed separable, and 
the bath is assumed stationary, the ground state evolu- 
tion during interval T can be neglected. Then 2D mea- 
surement essentially accesses the reduced evolution su- 
peroperator, Eq. and possibly also the more general 
superoperator 

W(^) (t, r) - trslW^-^' (t);7(3^) (r) W^eJ. (27) 

We will show below that this superoperator, together 
with the light properties, determines the way in which 
the molecule is excited in a general case, even at illumi- 
nation by natural light. 

IV. EXCITATION BY LIGHT 

In order to account for general light properties we will 
consider the problem fully quantum mechanically, and 
assume only deterministic evolution of the system wave- 
function. The Hamiltonian of the system reads 



+Hm^b+Hm~b+Hj^i-s+Hb-s+Hb-b+Hb.-s- (28) 

We have divided the system into a molecule {Hm), its 
environment or bath (Hb), the radiation {Hn) and the 
light emitting body (LEB) which produces it, e.g. Sun 
or laser medium (Hs). It seems reasonable to neglect 
a direct interaction between the molecule (together with 
its environment) and the molecules of the LEB. Conse- 
quently, the terms Hm-s ^-nd Hb-s can be disregarded. 
To make the treatment simpler we can also neglect the 
interaction between radiation and the molecular environ- 
ment, Hb-r- The assumption is that the energy of the 
molecular transition that is used to harvest light for pho- 
tosynthetic purposes is much larger than any of the tran- 
sitions in this environment and the two regions of the 
light spectrum can thus be treated separately. One can 
also assume that the part of radiation spectrum which 
would interact with the bath is simply filtered out, and 
the environment is kept at certain temperature by other 
means. 



A. Radiation entangled with the light emitting 
body 

An important special case is the one in which the ra- 
diation and the LEB is in equilibrium with each other so 
that the radiation is described by the canonical equilib- 
rium density matrix 

W^r^=J2 ^ z"" |A^A,)(A^A,|. (29) 

Here, \N\q) is the A^-photon state of the radiation mode 
with polarization vector ex and wave-vector q. As we 
have already noted above, the statistical concept of den- 
sity matrix will be replaced here with the concept of en- 
tangled states, so that we can describe the whole system 
by its state vector. Thus, we introduce a state vector 

I^W) = E E C^Aq(i)|A^Aq)|0^Aq(t)), (30) 
A<j 

in which the light is fully entangled with the states 
</'AfAq(i)) of the LEB . The LEB states have to fulfiU 
the condition 

{(t>N\qit)\(t)N'X'q'it)) ^ ^W'Sqq'SNx^N'^^, (31) 

SO that when the total density matrix of the LEB and ra- 
diation is averaged over the states of the body, we obtain 
Eq. (P^)) . W^'^'^ is recovered provided that 

|cArA,(<)P = " ■ (32) 

^Aq 

In the absence of the light absorbing body, the evolu- 
tion of the state |S(t)) is governed by the Hamiltonian 



H = Hm + Hb + Hr + Hs 



Hl^Hr + Hs + Hr^s, 



(33) 
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and 



UL{t) = e{t - to) 



xe^p^--iHs + HR + HR^s){t~to)} (34) 
is the corresponding evolution operator. 



B. Equation of motion 

For the subsequent treatment of the system dynam- 
ics, we introduce the interaction picture with respect to 
Hamiltonian operators Hm, Hb and H^, 

H^'Ht) = HM-sit) + HM-Rit), (35) 

where 

HM-Sit) = UlS)Ul{t)HM-BUB{t)UM{t), (36) 

and 

HM-B{t) = Uli{t)Ul{t)HM-BUL{t)UM{t). (37) 

Equation of motion for the total state vector in the in- 
teraction picture 

\¥'\t)) = ul,{t)Ul{t)Ul{m{t)) (38) 



thus reads 



||*<"W) 



- i + Hm-hH)) (39) 

n 

The solution of Eq. p9l) can be found formally as 



\^^'\t)) = exp+ j dT(HM-B{r) 



to 



H 



-h(t))}|*o). 



(40) 



We will assume that the system is initially in the state 
4'o) of Eq. dD). With this choice we have 



HM-B{t)\^o) =Q. 



(41) 



Further in this paper, we will assume weak interaction 
with the radiation, so that it can be described by lin- 
ear theory. Thus, we need to collect all terms in the 
expansion of Eq. (1401) which include one occurrence of 
HM-R{t). Thanks to Eq. (|¥T|) . however, all terms where 



HM-R{t) is not on the far right of the expression are 
equal to zero. Eq. therefore simplifies into a series 



|*W(t)) = |*o) - ^ j dTHM-R{r)\^o) 
to 



Z I 

^ JdT j dT'HM-B{r)HM-R{r')\^o) 



(42) 



Now we introduce a projector Ve that excludes the exci- 
tonic ground state \g) 



^e = ^|e„)(e 



(43) 



Applying this projector to Eq. (|^^ has only the effect 
of eliminating the first term of the series. Introducing 
abbreviations 



and 



E 

\S{t)) = -'- j dTHM-R{r)\^a) 

to 



|*«(t))=P,|*(^)(i)), 



(44) 



(45) 



we can write 



\^i'\t))^\S{t))-- I dTHM-B{r)\S{T)) 



to 



z r 

JdT j dr' HM-B{r)HM-B{r')\S{T')) + . . . . (46) 

to to 

It is possible to verify easily that this series is a solution 
of the equation 

= -'-^HM-B{mi'\t)) + ||5(<)), (47) 
with initial condition |^'i^''(io)) = 0. 



C. Pumping source term 

Eq. (|T7|) is an equation of motion for the excited states 
of an excitonic aggregate pumped by a source term 

\S'{t)) = ||5(t)) = -'-HM-R{m,). (48) 

Hamiltonian Hm-r will be assumed in the dipole ap- 
proximation, i.e. 



Hm-r w ^l ■ Erir), 



(49) 
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where /x is the transition dipole moment operator of the 
aggregate 



A* = ^d„|e„)(.g| + h.c 



(50) 



and Et is the operator of the (transversal) electric field 
of the radiation 

Erir) = -iJ2{ex-qf-q(.r)a{^ 

\q 



e\qfq{r)axq 



with 



(51) 



(52) 



Here, SI is a quantization volume. 

We consider a molecule much smaller than the wave- 
length of the light, so that e"' '^ is constant in the volume 
of the molecule. The origin of the coordinates can thus 
be conveniently put into the molecule yielding e"' '' « 1. 
The interaction Hamiltonian in Eq. (|^5)) then reads 

HM-Rit) = z^MAq(i)/A-q(0)at^(<) 

Xq 



«AiAq(t)/Aq(0)aAq(t), 



(53) 



where the creation and annihilation operators of the field 
are in the interaction picture with respect to Hamiltonian 
Hl, i.e. 



4qit)^Ulit)ai^ULit), 



axq{t) = Ul{t)axqUL{t). 



(54) 



(55) 



The transition dipole moment operator projected on the 
polarization vector of a mode \q appears in the interac- 
tion picture with respect to Hamiltonian 



f^\q{t) = Uli{t)n ■ exqUM{t)- 



(56) 



The evolution operator UL{t), Eq. can be rewritten 

as 



(57) 



where 



C/fl-s(t) - e(i - to)cxp+ I - f dTUl{T)ul{T) 



HR^sUR[T)Us{T)y 



(58) 



Since Hamiltonian Hs commutes with the radiation op- 
erators and 



Ul{t)axqUR{t) 



a-Xq, 



we have 



and 



axq{t) = aAq(i)e-^"'*. 
Here, we introduced slow oscillating envelops 

and 

axq{t) = U}j_g{t)axqUR-s{t). 



(59) 



(60) 



(61) 



(62) 



(63) 



Inserting these expressions into Eq. (|53p we can distin- 
guish two terms associated with the transition from the 
ground state |g) to an excited state \ea) with respective 
phase factors e*^""''""'')* and e*(""s+"'')*. While the first 
one will lead to a resonance excitation around ujq ~ uJag , 
the later term is oscillating fast and will therefore con- 
tribute very little compared to the former one. Thus we 
drop the fast oscillating part, and obtain the source term 
in the form 



Xq 



X axq{t)\g)\<S>B)\^a). 



(64) 



Using this form of the source term, we can find state into 
which the molecule is weakly driven by any type of light. 



D. Excited state dynamics under pumping 

So far we have treated the problem systematically us- 
ing the wavefunction approach. The time evolution of 
the system wavefunction is governed by Eq. (|39l) . To 
find the probabilities of creating population on and co- 
herence between certain excitonic levels \ea) we solve Eq. 
(|39ll formally. 



\¥JHt))= j drUM-B{t-r)\S'{r)). (65) 

to 

Here, we used the fact that \'^e{to)) = 0. Now let us 
evaluate matrix element Pab{t) ~ {'^ e{t)\'Pab\'^ e{t)) of a 
projector 



Vab = |ea)(e6|, 



(66) 
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which gives the probability of finding the molecule in 
state |ea) if a = 6, or characterizes the amount of coher- 
ence between states jca) and jcfc) if a ^ &. Note that we 
have removed the interaction picture, Eq. ([55)) . We have 

t t 

Pabit) ^^jdrjdr' J2 hi(fy<i'r 

to to ^I'^'l' 

x(So|alqWaA'q'(T')|So) 

x(eb|W(^)(i-T,r-r'K,,A'q'|ea), (67) 

where the evolution superoperator U'^'^^t ^ t,t — r') has 
been defined in Eq. (^T)) . 

In Eq. (|67p . the light is represented by a first order 
correlation function 

X (So|4g(T)avq'(T')|So) (68) 

(see e. g. Ref. [41]). which comprises all its relevant 
properties. We also denoted 

PXq.Vq' = ^MAq|5)(5lMA'q'- (69) 

The quantities Pab{t) are the matrix elements of the 
RDM {Pab{t) = {sb\p{t)\ea)) of the system which reads 

t t 

p{t)^ j dT j dT'U'^''\t-T,T -t') 

X E PAq,Vq'4q!vq(^,^')- (70) 
Aq, A'q' 

For a weakly driven system, Eq. (|70p has a very wide 
range of applicability. We will discuss its application to 
thermal light and pulsed coherent light in the following 
section. 



V. DISCUSSION 

Thorough discussion of excitation dynamics of molecu- 
lar systems excited by incoherent light was made in Ref. 
[stI . Molecular systems were considered without the bath 
effect which is however significant for light harvesting. 
Eq. ([70]) contains reduced evolution superoperator of the 
molecular system so that the state of the system created 
by the incident light depends on its reduced dynamics. It 
is not possible to consider a general case of such dynam- 
ics analytically, and we will therefore commit ourselves 
to some simple cases. 



In so called secular and Markov approximations (see 
e.g. Ref. [l^l) matrix elements of the evolution superop- 
erator governing the coherences take a very simple form. 
First, it is possible to separate the two time arguments 
in the superoperator W'-'^-' {t, r) so that 

U<^^\t,T)^U<^^\t0''3\T). (71) 

Since each coherence is independent of the population 
dynamics and of other coherences, the one-argument su- 
peroperator elements read 

^^il(0=e-— (^»+^^)*, (72) 

and 

iitgtit)=e-"^"^'-^-'. (73) 
Here the dephasing rate 

Ta = 7p + (74) 

comprises the pure dephasing rate 7^ and the rate Ka 
of depopulation, i.e. the sum of transition rates from 
state \ea) to other states. A simplified treatment of the 
populations is possible for the states that are only de- 
populated, i.e. no contributions to the population can 
be attributed to the transfer from other levels. They are 
found at the top of the energetic funnel of the antenna. 
For these states we have 

liitait) = e-^"*. (75) 

Eqs. (|72ll to ([75)1 neglect all coherence transfer effects, as 
well as possible coupling between the dynamics of popu- 
lation and coherence. 



A. Excitation of coherences by thermal Hght 

For an equilibrium thermal light the correlation func- 
tion -^AqA'q' (■'"'''"') depends only on the difference of the 
times r and r'. As discussed above, |So) represents the 
equilibrium of the system described by Hamiltonian Hl ■ 
The equilibrium density matrix is stationary, i.e. 

C/^(f)|So)(So|C/i(<) = |So)(So|, (76) 
so we can write 

4q,A'q(^-^') - |/Aqr'(So|4q(T-r')aAq(0)|So) 
Xe^"'(^^^')<5AA'<5qq- 

= |/Aqr'/Aq(r-r')5AA'V- (77) 
It can be shown that 

hq{~t) ^ il^it). (78) 
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Assuming some simple form of a light correlation func- 
tion, e.g. 



(79) 



we obtain for the populations 



r r 
^1 to to 



X e 



-r„(r-r')-^-«(<.^ag-c^„)(r-r') 



(80) 



Here, [p]ab = (calpleb). We utilized Eq. ((79)l and the fact 
that, by definition (see Eqs. (j68| and (|69|)). the time r 
corresponds to the action of the dipole moment operator 
from left, whereas time r' corresponds to the same action 
from the right. At long times t — to ^ oo this yields 



2 i^a + :f2^I\q[P\q]aa 

K'aiLOag-LOqr + iTa + ^J' 



E 



(81) 



However, neglecting the influence of environment as in 
Ref. [37|] yields 



Pant-to) = ^ 

\q 



(82) 



which grows linearly with time. 
For coherences we have 



1 



\q 



iUJab + (Fa + Tb) 



such a case fast enough to prevent averaging over many 
such spikes. Eqs. ([51]) and with large Ka describe 
just such a situation. The RDM created by incoherent 
light resembles in certain sense the one created by ultra- 
fast pulses; it represents a linear combination of excitons. 
The coherences in Eq. (155)) are however static at long 
times. 

In our demonstration we concentrated on a simple 
model assuming both Markov and secular approxima- 
tions to be valid. The presence or absence of coherences 
has no significance in such a case, and more involved the- 
ories of the RDM dynamics 0, have to be used to 
investigate the role of coherences in energy transfer pro- 
cesses by Eq. ([7n|. 



B. Coherent pulsed light 

In derivation of Eq. ([70)1 we assumed certain initial 
state I So) of the system composed of the light and its 
source. The condition that the light is in a stationary 
state, fully entangled with its source, has only been used 
to simplify the correlation function -f>q\/q/ (t, t') for the 
case of the thermal light. In a general case |So) will not 
represent an equilibrium state. It can indeed describe 
even systems such as a laser producing coherent Gaussian 
light pulses with some carrier frequency ujq and a width 
parameter A. If we in addition assume that the light is 
described by a single polarization, and that we consider 
the dynamics after one such pulse centered at time t — tq, 
the light is described as 



I'll . ,, (^"^of i^'-rg)^ , , 

(i) _ 7-„^ ^ ^34-) 



E 4',Vq'(T,r')=/oe 

Xq.X' q' 



i{uJag - LOq) + Ta + 



-i(uJbg 



(83) 



which turns into Eq. (|81|) for a — b (with additional as- 
sumption Ka = '2Ta). In case of no dephasing, the first 
fraction in Eq. ([55]) yields a delta function S^Ldat) W^ - 
Thus, for slow or non-existent relaxation due to inter- 
action with environment, the system is excited predomi- 
nantly into a state represented by a diagonal RDM, as all 
coherence terms are negligible compared to the linearly 
growing population. For fast relaxation, the coherences 
may be of the same order of magnitude as the popula- 
tions. 

The case of very fast relaxation is particularly interest- 
ing. It was suggested previously that coherent dynamics 
can be relevant for the in vivo case, because the fluctuat- 
ing light from the Sun corresponds to a train of ultrafast 
spikes ^JJ]. The relaxation of the antenna must be in 



Coherence element created by such light reads 
t t 

to to 



X e-^»(*-)-^^(*-')/oe-^-^pL, (85) 
where pJJ^ — -^{eb\iJ.\g){g\iJ,\ea). In the limit of ultrashort 

(t-tq)''' 

pulses when e ^ a6{T — tq) the pulse creates a 
pure state at tq, which then dephases as 

p,,(t) = e(i-ro)e-(r"+r^)(*-^°) 



(86) 



In case of a finite pulse and no dephasing our results 
coincides with those found in Ref. 13711. 
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C. Mediated excitation 

The major difference between excitation by the ther- 
mal hght and a coherence pulse is in the occurrence of a 
sudden event which populates a nearly pure state of the 
excited state band. Clearly, a single molecule interacting 
with an ideal continuum of radiation modes in equilib- 
rium does not experience such sudden events. Rather, its 
interaction with light corresponds to a continuous pump- 
ing, and the suddenness of the photon arrival is the conse- 
quence of our ability to register only classical outcomes. 
In order to register them we have to interact with the 
system and become entangled with it. Our experience is 
that macroscopic systems interacting with low intensity 
light can be used to detect single photons, and certain 
more or less definite times can be attributed to their ar- 
rivals. Interaction of a photon with a macroscopic de- 
tector yields a temporal localization of the arrival event. 
A mesoscopic system may play a role of such a detector 
(mediator) that provides its fluctuations to be harvested 
by dedicated nano sized antenna. Green photosynthetic 
bacteria, from which the photosynthetic complex FMO 
was isolated, collect light mainly by means of so-called 
chlorosoms (23 . |4^ . The chlorosom is a self-assembled 
aggregate of ~ 10^ BChls and carotenoids with very lit- 
tle protein. The typical dimensions the chlorosom are of 
the order of 100 nm ^|. It does not seem to be organized 
as an energy funnel [43, 44], and the energy transfer time 
between its main body and the base plate to which FMO 
complexes are attached is of the order of 120 ps [45)], i.e 
rather slow. The excitation in such a mesoscopic system 
may have enough time to become localized through in- 
teraction with the large number of the systems DOF and 
arrive at the FMO complex in a particle like, i.e. also 
temporally localized fashion. 

In this section, we will generalize our result, Eq. ([7n|. 
for when the excitation of the photosynthetic sys- 

tems occurs by transfer from another system. We will 
therefore assume that our molecule does not interact di- 
rectly with light, but is pumped in a similar fashion by 
another system. The source term, Eq. (j64l) . is then gen- 
eralized as 

\S\t)) ^ ^A{t)\g)\^B) 



Here, A = J2a.n \'^ri)\^g){^a\{g\ + h.c. is the molecule- 
mediator interaction Hamiltonian and the time depen- 
dence results from the interaction picture 

Ait) = ul{t)u\{t)AUA{t)UM{t). (88) 

We denoted the ground and excited states of the mediator 
by l^g) and |^„), respectively. The state of the molecule 



at long times is in analogy with Eq. ([70)1 
t t 

p{t)= [ dT [ dT'U^^Ht-T,T-T')Y,<ir)an'{T) 



X {UAiT)A{T')\^^,){MrMn'iT')). (89) 

The complicated two-point correlation function in Eq. 
(j89p results from the pumping of the mediator similarly 
to the direct pumping of the molecule in Eq. ((70)) . A 
mesoscopic system especially when excited will, however, 
always exhibit fluctuations which will prevent the corre- 
lation function from having a simple smooth dependence 
without recurrences. Such recurrences can temporally lo- 
calize the excitation events of the molecule. In such an 
excitation regime, when coherent dynamics from differ- 
ent excitation times do not interfere, optimization of the 
FMO's energy channeling capability for case of initially 
coherent states would be an advantage. 



D. Outlook 

More research into specific forms of both the light cor- 
relation function for different situation that may occur in 
vivo, and the analogical interaction of systems like FMO 
with mesoscopic antennae is clearly needed. Ultrafast 
spectroscopic experiments play a pivotal role in this re- 
search by yielding information about the system response 
to the light. To conclude on the utility of coherent dy- 
namics for the function of the photosynthetic system is, 
however, only possible by taking into account the prop- 
erties of light at the natural conditions, for which the re- 
sults of this paper provide means. If the coherent dynam- 
ics observed in some photosynthetic chromophore-protein 
complexes has a significance for their light-harvesting ef- 
ficiency, and these systems evolved to optimize it for their 
corresponding ecological situation, it can be expected 
that the properties of at least some parts of the pho- 
tosynthetic machinery would be tuned to the fluctuation 
properties of their source of excitation. For plants and 
some bacteria this may be the Sun light, others like FMO 
complexes could be expected to be tuned to the proper- 
ties of their associated chlorosoms. 



VI. CONCLUSIONS 

In this paper we have discussed dynamics of a molecu- 
lar system subject to external pumping by a light source. 
In particular we have considered excitation by thermal 
light, by coherent pulsed light and an excitation through 
a mesoscopic antenna. With a completely quantum me- 
chanical treatment, we have derived a general formula 
which enables us to study the effect of different light 
properties on the photo-induced dynamics of a molecular 
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systems. This formula naturally contains the system 
environment interaction contribution to the excitation 
process which enters via appearance of the reduced den- 
sity matrix dynamics. We show that once the properties 
of light are known in terms of a certain two-point correla- 
tion function, the only information needed to reconstruct 
the systems dynamics is the reduced evolution superop- 
erator, which is in principle accessible through ultrafast 
non-linear spectroscopy. This conclusion applies to any 
type of light and makes thus the results of ultrafast spec- 
troscopic experiments universally relevant. Considering 
a direct excitation of a small molecular antenna we found 
that excitation of coherences is possible due to overlap of 
homogeneous line shapes associated with different exci- 
tonic states. These coherences are however static and 
correspond to a change of the preferred basis set into 
which the system relaxes from the one defined by the 
bath only, to the one defined by the action of both the 
light and the bath. When an excitation of a photosyn- 
thetic complex mediated by a larger, possibly mesoscopic, 
system is considered, the complex can harvest fluctua- 



tions originating from the non-equilibrium state of the 
mediator. Fluctuations of the mesoscopic system such as 
chlorosoms may time localize excitation events of the en- 
ergy channeling complex, and to excite adjacent energy 
channeling complex coherently. It is likely that in such 
a case the properties of energy channeling complexes like 
the well-know Fenna-Mathews-Olson complex would be 
specially tuned to the fluctuation properties of their as- 
sociated chlorosoms. 
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